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Abstract 

Hyp* ("t roi nih tit: s are a gent rahzation to second ordt r (r n sor Ji elds <>) fh< . :>nrf n- 
tionai >tn a mimes used in vector field visualization. .Is nppustd to point iron s common!'/ 
it** d in rt*U'ihzin >j tensor fields, h ijpcrstrcamlines furm a > >ntm nou> n pr> s* ntation -f 
th f ooniplt tf ft nsor information along a three-ihmension ol path. This technique is ns< ful 
in n^nahzing both > i/mmtfnc and nnsymmetru: thn e-dune n<i<>nal tensor data. .A » * ral 
>rawplts of tensor field visualization in solid materials and thud flows arc pmnded. 


keywords - scientific visualization, multivariate data visualization, tensor field > i~ 
sualization, flow visualization. 



OWQJNAi PAG* m 
QUALfTY 


1 





1 Introduction 


\ I . my pay-. -ml :>:• ■ ! i*;i rr**itT !y vi-uali/ed incompie'.-iy in *. rm> of v- c*, • >r 

1 1 u i .1 ’ ■ < • 1 l -i , a •’ X 1 '« t | T * * M [ SOT I ! I S f ) 1 <1 V t * 1 1" ! i ! 1 1 ( 1 1 1 f ' S . I * [ ,f >V i < 1 ■ ! : a ' < ' , c X ’ i - f . 

hod< >i» igy fi ir * ; - mdi / i n g 3 - 1 ) - e r< > n d order T em-or dat a ran mad Jn new ] u -;u h f - ' . . > - 

problems. Second order fen>or fields are fundamental fn many disciplines of’iM'uiic^n'id 
physical -chuices. St resses and -train- in soiids. for I'xam pie. an* tensor fends. In : ! i : i ■ i Aom. 

reuses, viscous stresses. rate-of -a rain, and momentum ■ ran-fer< are all described in T m 
of tensor data. Also, t ho steady-state Xavier-Sf okes *m | ua t ions describing gas flows invAp.m 
oiil y out? quantity: a tensor field railed moment um Jinx tl< nsity. Fable 1 gives the < i efi ji i ' : • o : 
of some common 3-1) tensor fields. It shows that tensor data are very rich in information 
content: they include diverse physical quantities -uch as pressure, kinetic energy den-iyy. 
mass density, velocity, and derivatives of the velocity held. Visualizing tensor fields proved-- 
.l correlation between these quantities. 

As a result of this wealth of highly multivariate information, tensor vmuaii/at j» m i- a 
ehallenging enterprise. Indeed, a 3-1) second order ?en-nr field T consists of a 3x3 array *u 
.-(‘alar functions {7^}. i.k = l . 2. 3 defined over a 3- 1) domain 1 . Independent visualization , .f 
these nine functions is possible but meaningless. Fhis article presents a methodology ba-ed 
on the concept of a hyper streamline* which is the simplest continuous tensor structure th.it 
can be extracted from a tensor field (as opposed to many other scalar or n (tor feat tires >. 

In Mm i ; e \ t section, hyperst reamlinos are introduced for the particular case of .-// m~ 
no trie tensor fields U = { f A } whose individual c< an ponents are related to each other by 
I lor i*k — 1.2,3. As seen in Fable 1. symmetric tensor fields are very common in 

fluid how studies. Then, a structural depiction of symmetric tensor fields is derived from 
the representation of a large number of hyperst pm mii tm-. Finally, a methodology to visu- 
alize uusynwH trie tensor data is provided by encoding' an additional vector field along the 
trajectory of the hyperstreamlines. 

‘ t s.iful in tXrmation about mn>or fields can lx* found in a- i 1 . 
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I able 1: I'-Misor fields in fluid !lows. p — pressure. p = mass density, r, and - v> ■ 
component'. and .7 — viscosity. is the Kronecker ye m I10I. nn - u n -yinm* eric. - 
symme» rir . 


2 Symmetric tensor fields and hyperstreamlines 

Symmetric tensor fields U = {fYv} — {f\t} aro very common in physics and : 1 ll i t : * 1 ► 1 r; r • -2. 
and their visualization is an important problem per so. Also, it is a necessary - : ■ ■ p fn 
display of general unsymmetric data, winch will be addressed in the last -ecAon *. -!’ r M- 
art icie. 

A symmetric tensor U may be described by three orthogonal vector f I o ] d - . Inn* * I 
1 1 0 m at every point T three real eigenvalues A'W i — 1 , *J. or A ordered according to 

A (l) > A iJ > > A n > 


as well as three real and orthogonal unit eigenvectors W- f 1 ] . We consider the three orthog- 
onal vectors'll jriven by 

- (l) = A (t) r (0 . 1 . 

IJoca u s*' of the particular ordering of the eigenvalues, we refer to T' {[ - as the major < '<;• ; - 
( < 'nr. 7 l ^ as the medium eigen vector, and as 1 lie minor eigenvector, 

\ Uuaiiziug U is fully equivalent to visualizing simultaneously the three vector field- 7 ' 7 
since they include all the amplitude information (the eigenvalues A (: M and all the din rt : ura! 
information (the unit eigenvectors “bl) represented in matrix notation by the component- 
Furthermore, visualizing the three vectors 7 !:! allows one to understand the b<’!iav:«-r 
of the six indepeiulont components Ij^ with little or no training [ 2 ). 
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! l . i v i : i g ‘i • r ’ .1 r;:,''; ;>.i i . i \ • •- * »e t 1 1 r* v* t u *r - 1 i > r I ! , 1 1 > ► • r - ** ■ >< > r p - ■ > { , . i V , ; ; , * . ■ . . • ^ 

omp'-m. 1 i,. mb > rr £! a r :* >n ;i! a u: vmi * * x ’ . 1 1 h >u :>ut ' : i * • : r • ; 1 - < ' » 1 ■ . .* 

d Mm underlying coioitiuiry « »f Mm data : i i * i i L l'- >r t his ma-om I ’) ; * * k ’ : l - * ns -V a 1 -■;- 
jit Id hn* s which are m reamliims of « urn of ‘ lit' eigenvector • i < ' i < J s ■ *. Whim m:i| urn.-;/: :.g • :.r ,j 
rout inui tv. tensor t i ♦ *i t i lines are only nrior nans and represent nnlv parMaliv ’ lm mu -<«r 
fiold (this is analogous to representing a vector field by examining oiilv mm of jts ;i,ar 
component s i. 

What in more desirable is a lint ftvsor >rnn that r*'presen f - all the oei-or in form.i ■ «n 

along a T 1) pat h in space or. ♦'<p.ii vahnitly. that encoder a <*■ >nt inuous di- f r; but ion of ip-« -id- 

along a given t ra ject< mv. For tli’m purpose v,v generalize Mm viator notion of a -tr-uimidm 

to t ho tensor concept ot a hfjp* r>/;v timiint : 

a geometric primitive of finite size sweeps along one of the eigenvector imfd- 
T 7 '^ while >t retching in the transverse plane under Mm combined action of Mm 
two other orthogonal eigenvector fields. h!ie .surface obtained by linking the 
stretched primitives at the different points along the trajectory is called a 
;>f rstn amlint and is color coded by means of a u-er-dej; n^d function of the three 
eigenvalues, generally the amplitude of the longitudinal eigenvalue. 

The color and trajectory of a hyperst reamline fully represent the longitudinal eigenvector 
field and the cross-section encodes the two remaining transverse eigenvector fields. Thus. 
Iiyper-treamlines form a continuous representation of the whole tensor data along the trajec- 
tory. Ilyperstreamlinos are called major , medium, or whtor depending on the longitudinal 
eigenvector field “td that, defines the trajectory. 

figures 1 and 2 illustrate the properties of hy p»u>r reamlines for two elastic -tress fields 
in a semi-i nii nite steel solid, and are described in detail below. The color scales for these 
and every other figure in this article are shown in 1 Mg. 2. 
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Figure 1: Stress tensor induced by two compressive forces. Hy perstreamlino t ra p ' r t oil* s 
(top): minor t u bos . medium and major helices D xjttunij ; solenoidal property of a minor 
tube (right). Color scale A of Fig. 2 is used. 


2.1 Trajectory of a hyperstreamline 

Hy perst reamline trajectories correspond to Dickinson's tensor field lines o : these patterns 
of lines -how, for example, how forces propagate in a stress tensor held, and how the 
momentum is transferred in a momentum flux density tensor field. Figure l(top) illustrates 
tins phenomenon in an elastic stress tensor field induced by two compressive forces on 
the top -urface of the cube. The lines propagating upward are along the most compressive 
direction . the minor eigenvector 7;^ ). and converge towards the regions of high stress where 
the [Drees are applied. Note the sudden divergence of dose trajectories on each side of the 
plane >,[ -y mmet ry. Similarly, trajectories along the two other eigenvectors delineate a 
surlace -tiown near the horrom lace ot tlu? cube. Lids surface i> everywhere perpendicular 
to the mo-r compressive direction. 
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2.2 Cross-Section of a hyperst reamline 

IF. :•< : ; : , 1 1 : . * • ^ ;i a- : u r ' : m r characterized by t !m m’omet ry <>l \ ! i * m r cr< mm^'cMon. i m. ’ ;m 

jm< mm* primi* iv*> ‘a. it sweeps along • he *ra jam < »ry. Wa o ui>id»T *wo yvpes rtf p r i : : i d 1 v*~; 
i a eov/t mm ch»‘s into mi < ■ 1 1 i p ^ * 1 wloh 1 . t nil t li.it mummim a a vpam* r*u;:;- 

i.:m caimd ,* /aa*; am! 2» a ao>.>> * ii.tr ganarates a 1 1 y [)♦.* r s t reamlina called a In Us. I m ; : • ' 
F bottom) allows t wo minor tubes propagat ing upward as well as four madium ami : : : o b>r 
Indices in Mm stress i< i n>or field corresponding to Fig. It top ). In a tube, the principal ax- ^ 
of *’ach elliptical rross-sertion are along the transverse eigenvectors and have a length pro- 


portional ’o the magnitude of Mm transverse eigenvalues. The same property holds for 
helix, whose arms are proportional to the transverst' ajgen vectors i helices owe their name 
to the spiraling pattern of their arms that can be observed in some cases). In this manner 
both directional and amplitude information are encoded along the trajectory. The local 
sign of the transverse eigenvalues can be detected by examining the singularities in Mm 
cross-sect ion of the hyperstreamline. Indeed, the cross-section reduces to a single line or a 
point wherever one of the transverse eigenvalues changes sign. 

1 ubos and helices encode the same information about the tensor field, but some a>pecm 
ol i he data are bettor perceived with one hyperstreamline than with the other. Tubes, for 
mmance. show better where t lie tensor is degenerate in the transverse plane, since recog- 
nizing that an ellipse is circular is easier than comparing the length of two perpendicular 
line segments. Further, if the tensor field is transversely degenerate in a whole region of 
space. Indices are not adequate since in this case the direction of the transverse eigenvectors 
i." not determined. Helices, on the other hand, provide better cues for perceiving preciodv 
Mm direcMons of the transverse eigenvectors. 

f our different stages of a minor tube in a stress tensor held are displayed in Fig. 2. Tim 
teimor field is similar to that of Fig. 1 but an additional tension force is added. In the top- 
lot t . the cross-section is circular, and the transverse stresses are equal in magnitude. The 
top-right shows an increasing anisotropy of the transverse stresses together with a decrease 
of the longitudinal eigenvalue (color). In the hottom-hdd. Mm cross-section is reduced to 
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Figure 2: Four different stages of a minor tube in an elastic stress tensor field. ('olur -c.ih 
A is ii soil. 

a straight lino: one transverse eigenvalue is zero and the stresses are locally 2- IF In A.* 
bottom-right. the stresses are 3-D once again: the eigenvectors undergo a rapid ru\ ; A<e 
and a -ubstauMal stretching which reveals an important gradient of Aiear and pre--* a re •• 

‘ 1 i * 1 Y< l ga. Hi . 


Degenerate and singular points. Computing hyperstreamlines is complicated because 
degeneracies can occur along the trajectory at and in between the sampling points requested 
by the adaptive integration algorithm. We assume that the tensor field is smooth, i.e. Aiat 
the direction c>f the longitudinal eigenvector is not likely to vary by more than a us»t- 
pr-'denin'd angle between two successive sampling points, unless the trajectory just civ.— A 
a degeneracy involving the longitudinal eigenvalue. In tins case, we search for t he degeneracy 
bei ween the hist two sampling points and, if found, terminate the curve there. Wo can then 
jump t he degeneracy and continue integrating in a selected eigen direct ion. The points where 
the transverse eigenvalues vanish are also detected and included to the curve in order not 
to miss a singularity of the cross-section of the hypei>t reamline. 
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SoltMioiilal tensor fields 
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Fxampms of snmnoidal vector fields include the vorticitv 


nr ’lie vmoci 


compressible Hows. Much of the structure of sohumidul vector fields <. 


'ail ! > e * 1 x T I [ ; : ] t : * 


their property of having a constant fltix inside a st ream? 'fee ’ ] : the magnit ude <*f * !,»• v. ,r 
field must increase locally in regions where the streamlines converge towards each .md 

decrease where the streamlines diverge from each other. 

lly analogy, we define a tensor field U as >o!f noidul if it satisfies 

iJx ‘ 

for k — l.’i.d which implies that the three vector fields obtained by multiplying U by anv 
three constant orthogonal directions are solenoidal. 

Sohuioidal tensor fields are not rare mathematical objects. They are fundamental to 
fluid and solid-state mechanics. For example the stress tensor rr, k in solids at. rest An 
regions where no external forces are applied) and the momentum flux density tensor IF.. 


in gravity- free steady-state fluid flows both satisfy Kqn. 2 and are solenoidal tense 


tsor n<uu- 


umpt.ion of no gravity is common practice v. hen computing gas flows). Ihiwevr. 


i . * ‘ Mre>s and viscous stress tensors in fluid flows a 


in >t m > mnoi dal . 


Ilyper-t reamlines of solenoidal tensor fields have a c< uivergonce/di vergence propertv 
analogous to the property of streamlines in solenoidal v.vh.r fields. More precisely. jf .V 1 A 
the longitudinal eigenvalue of a hyperstreamline along Am eigenvector T' n and if A'- ' ami 
\ are the two transverse eigenvalues, Fqn. - is <mmvnhm! to 


A 1 *) - Ah ,( A t:| - \ {J} ) - A' 
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figure d: Rake of major mhos of tho inoinon t n m :!ux density fonsor in tho how : > a r_ an 
ogive cylinder, (\jlor -('ale A of Fig. 2 is u-ed. 

where A 1 ' 1 ' is the derivative ul t. he longitudinal » ■! gen va 1 ue a 'It/ruf fin fi'njt r/ory i ,i . /v 
A, ; . are geometric lactor- that are positive 1 it neigh m u’ing hvperst roannines ahum ’c <: ' 1 cnn- 
verge in the corresponding transverse direction.-, j or e and negative otherwise. Thus, when 
the longitudinal eigenvalue of a major hyperst reaimiue increases during its propagation 
i A m > O n neighboring major hyperst reamline< converge inwards each other. Anv diver- 
gence m one * rau^veroM‘igeii(iirection must he compensated by a stronger convergence in tie* 
other t «;a n.sverse (‘igendirecrion. The opposite property hold- true for minor hvperst ream- 
lines: an increasing longitudinal ( minor i eigenvalue i-urnAiNs with a global divergence of 
neighboring minor liy per.-t reamlines. Iftiiev converge in mi»' * ran-ver-e eigendirect ion. thev 
must diverge more strongly in the other transverse - ig-umlo'e: ;on. ( 'onversely, a decreasing 
longitudinal ‘Ugcuiv;tl ue corresponds to diverging ';:a;nr in. ;wr>t reamlines and converging 
minor hyper.-: reamlines. 

It tollows. t lien . that encoding the longitudinal ejumivalue into the color of a hvper- 
st reamline in a sobuioidal tensor held gives inlnrmat n ui aboijt the behavior ol neigh b< >ring 
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hvpmsrr< i amlinos. The figure at t ho beginning of this article illustrates this property fur 
major tubes of 'he momentum flux density Tensor in the flow past a hemisphere cylinder, 
hhe decrease of the longitudinal eigenvalue (color) is accompanied by diverging trajectories. 

Nie former property explains why the minor hyperstreamlines in Fig. 1 converge to- 
wards the applied forces (quickly decreasing longitudinal eigenvalue) and why the major 
hyperstreamlines propagate mostly parallel to each other with an almost, constant color. 
A close view of the sudden divergence of minor hyperstreamlines on one side of the plane 
ot symmetry is given in Fig. 1( right). The local divergence of minor trajectories creates a 
sudden increase of the longitudinal eigenvalue counterintuitive to the notion that the minor 
eigenvalue should decrease uniformly when approaching one of the two applied compressive 
forces. 

Figure 3 shows a rake of major tubes of the momentum flux density tensor in the flow 
past an ogive cylinder. The air flow comes in from a direction 5° to the left of the ogive axis 
and vortices are created in the wake of the body. Major tubes that become entangled in the 
vortices undergo a fast decrease in color while diverging from each other. In other regions 
of the flow, the color is constant and in some places it even increases slightly from orange 
to red. In these regions the apparent divergence of the tubes in the direction parallel to the 
surface of the body is compensated by a stronger convergence in the perpendicular direction. 
Both divergence and convergence exactly compensate each other in the tail between the two 
vortices. 

2.4 The reversible momentum flux density tensor 

A specific example of fluid flow analysis illustrates how hyperstreamlines may be used to 
correlate several different physical quantities. For the reversible part of the momentum 
flux density tensor, IT^ (see Table 1), one may correlate pressure velocity direction 1,., 
and kinetic energy density k [2]. Indeed, the major eigenvalue of is A^ = p + 2k and 
the corresponding unit eigenvector is the velocity direction I L ,. The other eigenvalues are 
degenerate (A^ 2) = A^ — p) in the whole space. It follows that only major tubes can be 
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Figure l: Reversible momentum flux density tensor in the flow past a hemisphere cylinder. 


Color scale A of Fig. 2 is used. 


used. Their trajectory is everywhere tangent to the velocity direction l v and their cross- 
section is circular, with a diameter proportional to the pressure p. The color of the tubes 
is determined by the function 


color ^ 


Ad) _o.5(A (2 > + A< 3 >) 


(4) 


which represents the kinetic energy density k. Thus, the trajectory, diameter, and color of 
the major tubes encode the velocity direction, pressure and kinetic energy density, respec- 
tively. 

Figure 4 shows \V tk in the flow past a hemisphere cylinder. The direction of the incoming 
flow is 5° to the left of the hemisphere axis. The detachment at the end of the cylinder is 
clearly visible. The pattern of hyperstreamlines indicates that the momentum is transferred 
from the tip of the body to the end fairly uniformly with a globally decreasing kinetic energy 
as shown by color variations. However, there is a sudden change of kinetic energy (color) 
and pressure (diameter) associated with a significant variation of the direction of the first 
five tubes. 
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2.5 Color coding schemes 

Usually, color encodes the longitudinal eigenvalue in order to represent the whole nei>or 
data along the trajectory. In practice, the color coding scheme can he modified to reveal 
other aspects of the data. An example is the color coding function of Fqn. 4 which allows 
decoupling of pressure and kinetic energy density when visualizing the reversible momentum 
transfers in a flow. For stresses in solids or viscous stresses in fluids, color can he usod to 
discriminate between compressive and tensile directions in the cross-section of a tube. This 
is done by coloring the tube according to 

color ~ cos ( ^) i o ) 

where gr is the angle between the normal n to the elliptical cross-section of a tube and 
the force / = U7T acting on it. Figure 5 represents the same minor tube as in Fig. 2. but 
colored according to Eqn. 5. Red corresponds to y = 0° and indicates that the corresponding 
directions 71 are in pure tension. Blue indicates purely compressive directions (y — 1M)°«. 
and green reveals pure shear = 90°). 

When using the color function of Eqn. 5 for other tensor data, the meaning of compres- 
sive and tensile directions is lost. However, this scheme encodes the sign of the transverse 
eigenvalues: a principal direction of the elliptical cross-section is red if the corresponding 
eigenvalue is positive, and blue otherwise. 

3 Structural depiction of symmetric tensor fields 

Two factors limit the practicality of hyperstreamlines: 1) the resulting display depends on 
the initial conditions of integration and 2) a large number of hyperstreamlines produces 
visual clutter. The same problems arise in 3-D scalar and vector field visualization. For 
example, when visualizing a scalar field with isosurfaces, the final image depends on the 
particular isosurfaces chosen and only a few of them can be displayed simultaneously. Also, 
the conventional streamlines used in vector field visualization lead to a display dependent 
on the initial conditions of integration and the presence of too many streamlines clutters 
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Figure 5: The minor tube of Fig. 2 colored as a function of the normal force. Color scale B 
of Fig. 2 is used. 

the image. In the latter case, these problems are overcome by algorithms that extract 
automatically the vector field topology [7, X], These algorithms can be seen as a way of 
coding the collective behavior of a large set of vector streamlines. Analogous to these vector 
techniques, a structural depiction of tensor data can be obtained by coding the collective 
behavior of a large number of hyperstreamlines. 

Consider the collection {H S^} of hyperstreamlines propagating along the eigenvector 
field as given by Eqn. 1. Important features exist in both the trajectory and the cross- 
section of these hyperstreamlines. For example, the locus 

A (,) = 0 

is the set of the critical points 2 in the trajectory the hyperstreamlines {//S* 1 *}. Further, 
the surface 

A (j) A (fc) = 0 

2 At a critical point of a vector field, the magnitude vanishes ,iml the direction of the streamline is locally 
undefined. See Ref. [9] for a complete discussion of this topic. 
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Figure 6: Structural depiction of t he stress tensor of Fig. 1. 


whore A^'l and A {v > are the transverse eigenvalues, is the locus of points where the cross- 
section of the hyperstreamlines {/f5^} is singular, i.e. is reduced to a straight line or a 
point. In general, a surface of constant eccentricity is the locus of points where the cross- 
section of each hyperstreamline in {HS^} has the same shape, regardless of its orientation 
and scaling. In particular, the locus 

± A<*> = 0 

is the set of points where the cross-section degenerates into a circle (zero eccentricity) 3 . 

A structural depiction of the stress tensor of Fig. 1 is given in Fig. 6. The yellow surface 
is the locus of critical points of the medium eigenvector T^ 2 ) and the green surface represents 
the critical points of the major eigenvector T 7 * 1 ^. On both of these surfaces, the cross-section 
of each minor tube (four of them are shown) reduces to a straight line. On the blue surface, 
the transverse eigenvalues are opposite to each other and the cross-section is circular. Below 
the yellow surface, both transverse eigenvalues are positive and every transverse direction 
in the cross-section of the minor tubes is in tension. Above the yellow surface, the medium 
"The locus X J} + A' fe ' — 0 was omitted in Ref. [2] . 
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-’nmnvnlne becomes negative and ^ome transverse directions are in tension while others are 
c( impressive. Inside the green surface, however, every transverse direction is in compression. 

3.1 Stress and viscous stress tensors in fluid flows 

Another example of .structural depiction is given for the stress tensor rr lk and the viscous 
stress tensor a[ k in fluid flows. As shown in Table 1, these two tensors differ only by an 
isotropic pressure component, implying that the unit eigenvectors of both tensor fields are 
identical. However, the eigenvalues of cr x k are equal to the eigenvalues of cr' lk minus a large 
pressure component. Table I also shows that visualizing a[ k is equivalent to visualizing the 
rate-of-strain tensor in incompressible flows. 

Hyperstreamlines of the stress tensor in the flow past a hemisphere cylinder are shown in 
Fig. 7 ( top ) (the flow is the same as in Fig. 4). The major tubes in front are along the least 
compressive direction Their trajectory shows how forces propagate from the region in 
front of the cylinder to the surface of the body. The cross-section of the tubes is circular, 
indicating that the pressure component of the stresses is dominant, as expected. However, 
t he viscous stresses close to the body create a slightly anisotropic cross-section. On the 
yellow surface, the eccentricity is equal to 10%. 

The helices are along the medium eigenvector field. They propagate mainly parallel to 
the cylinder surface and the orientations of their arms indicate a fairly constant direction 
of the two transverse eigenvectors. The third helix exhibits a more complex behavior, 
suggesting that the stress tensor is less uniform in the region of contact between the tubes 
and the body than in other parts of the flow. 

Figure 7( bottom) shows the viscous stress tensor a[ k in the same flow. As expected, the 
trajectories are similar to those in Fig. 7(top), but removing the large isotropic pressure 
contribution dramatically enhances the anisotropy of the cross-section of the tubes. The 
surface corresponds to a constant eccentricity of 90% and is crossed twice bv each tube. 
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Figure 7: Stress tensor (top) and viscous stress tensor (bottom) in the flow past a hemisphere 
cylinder. Color scale A of Fig. 2 is used. 

4 Unsymmetric tensor fields 

fly perst reamlines an' useful in visualizing symmetric tensor fields whose three eigenvector 
holds given by Eqn. 1 are real and orthogonal. However, visualizing unsymmetric 
3x3 tensor fields T = {Tj*.} is more difficult because their eigenvectors are generallv 
complex and not orthogonal. A 3-D vector field visualization technique, the stream polygon 
[10], can reveal aspects of vector field gradients. It does not apply, however, to other kinds 
of u asymmetric data, and the tensor information, even for vector gradients, is only partially 
rendered. 

In this section, we show that it is always possible to decompose unsymmetric tensor data 
into two components: a symmetric tensor field and a vector field. The symmetric tensor field 
is visualized with hyperstreamlines as before. However, in order to represent the complete 
(unsymmetric) tensor information, we need to encode the additional vector field along the 
trajectory. Depending of the physics involved, two reductions of the unsymmetric data are 
possible: a symmetric / antisymmetric decomposition or a polar decomposition. 



4.1 Symmetric / antisymmetric decomposition 


The tensor field can be decomposed into the sum of symmetric and antisymmetric cumpo- 
nent.s according to 

T 4 - T* T - T* 

T = 

2 2 

where T' is the transpose of T. The antisymmetric tensor has only three independent 
components that form a vector known as the axial vector [1. 11]. For instance, the velocity 
gradient in fluids is the sum of the rate-of-strain tensor c x 'k (symmetric) and the rate-of- 
rotation tensor (antisymmetric) which is half the vorticity vector. 

Figure O(top) shows aline tensor icon for unsymrnetric data based on this decomposition. 
A hyperstreamline is integrated along one eigenvector field of the symmetric tensor 
component and is color coded either according to the longitudinal eigenvalue or as in F<pi. 
5. An additional ribbon is added outside of the tube surface in order to represent tin' 
axial vector. The ribbon position and width encode locally the vector component which is 
perpendicular to the trajectory. The color of the ribbon maps the angle between the axial 
vector and the direction of propagation of the tube according to color scale B of Fig. 2 ‘ red 
is parallel, green is perpendicular, and blue is antiparallel). In Fig. 9(top) for example, 
color shows that the vector field is everywhere close to alignment with the direction of 
propagation. It is, however, not exactly aligned since the ribbon has a finite width. 

When visualizing the velocity gradient in fluid flows, this icon shows the position of the 
vorticity with respect to the principal strains, which is an important factor for understanding 
turbulence [12]. 

4.2 Polar decomposition 

An alternative reduction of the unsymrnetric data is the polar decomposition [11], which is 
a generalization to tensors (or matrices in general) of the usual decomposition of a complex 
number into the product of an amplitude and a phase. Assume as in Fig. 8 that the tensor 
T at a given point x maps the vertices of a cube from an initial state to a final deformed 
state. This global deformation can be decomposed into more elementary transformations. 
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Figure 8: Polar Decomposition of u asymmetric data. 

For example, one can first stretch the cube by a tensor U and then rotate the stretched 
rhomboid by an isometric transformation Q in order to reach the final state. Alternatively, 
one can first rotate the cube and then stretch it by the tensor V. 

Mathematically, these are two equivalent ways of decomposing the unsymmetric data 
into the product ot a stretch tensor U or V (the amplitudes) with an isometric transforma- 
tion Q (the phase): 

T = QU = VQ (<>) 

where both U = VTT and V = \/TT t are symmetric positive definite tensors, i.e. sym- 
metric tensors having real and positive eigenvalues, and Q = TU' 1 is an orthogonal tensor. 
It. can be shown [11] that this decomposition is unique wherever detT ^ 0, i.e. there is a 
one-to-one correspondence between the matrix T and the set of matrices {Q, U, V}. We 
will explain below how we handle points where det T = I). From now on, we restrict our 
discussion to the first decomposition in Eqn. 6 without loss of generality. 

Figure 9( bottom) shows a line tensor icon for T based on this decomposition. Tim 
symmetric tensor U is represented by a tube along one of its eigenvectors. In regions when' 
detT > 0, the isometric tranformation Q is simply a rotation and is characterized by an 
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Figure 9: Examples of line tensor icons for two different unsy unmet ric tensor fields: symmet - 
ric / antisymmetric decomposition (top) and polar decomposition (bottom). Color scales A 
and B of Fig. 2 are used for tubes and ribbons respectively. 

angle 6 (0 < 8 < x) and a unit axis of rotation 77. Thus. Q can be represented bv the vector 

7 (l > =*77 

However, a rotation of angle 8 about the axis 77 is physically identical to a rotation of angle 
2tt - 8 about the axis —77. Thus, Q is also equivalent to the vector 

“1“) = ( 2 t — 8 )( —77) 

To visualize Q, two ribbons are added that represent the vectors 7^ 1J and respectively. 
Note that in regions where det T < 0. the transformation Q involves an additional inversion 
in the direction of the rotation axis. . v row of white pearls across the ribbons marks the onset 
of the inversion during the propagation, and a row of black pearls indicates its cancellation. 

Singular points. In some points of the trajectory, the vectors 7 (1 ) and 7^ 2 ) may not 
be defined. We then simply interpolate them between adjacent points in order to avoid 



discontinuities in the ribbons. These singular [joints occur a) at the pearls where dot U = 
det T = t) i U is defined but not invertible and Q can not be computed) and b) where 
Q reduces to plus or minus the identity matrix (the rotation axis 77 is undefined and T is 
locally sy u 1 me i r i c with all eigenvalues having; the same sign). Thus, the assumption is that 
-insularities are isolated points along; the trajectories. This approach fails only if there is 
an entire sub volume where condition a ) or b) occur. In the latter case, however, simple 
hyperstreamlines are useful since the data is symmetric. 

5 Conclusions 

The wealth of information contained in second order tensor data is extracted and rendered 
as hyperstreanilines. By representing continuously both the amplitude and the directional 
information typical of tensor data, hyperstreamlines reveal much of the physics involved 
in complicated processes that are otherwise only partially visualized in terms of vector 
or scalar functions. Hyperstreamlines are the simplest continuous tensor structures that 
can be extracted from symmetric or unsymmetric tensor fields, and coding their collective 
behavior is a first step in obtaining a structural depiction of tensor fields analogous to 
extracting vector field topology. Future work must be carried out to obtain more advanced 
structural depictions. To this aim. it might be necessary to focus on specific tensors each at 
a time and to use the known underlying physics and the resulting tensor properties within 
the framework of this article. 
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